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ABSTRACT 



M. J. Holman 2 , M. J. Lehner 4,1,2 



We conducted a search for occultations of bright stars by Kuiper Belt Objects (KBOs) to estimate 
the density of sub-km KBOs in the sky. We report here the first results of this occultation survey of 
the outer solar system conducted in June 2007 and June/ July 2008 at the MMT Observatory using 
Megacam, the large MMT optical imager. We used Megacam in a novel shutterless continuous-readout 
mode to achieve high precision photometry at 200 Hz. We present an analysis of 220 star hours at 
signal-to-noise ratio of 25 or greater. The survey efficiency is greater than 10% for occultations by 
KBOs of diameter d > 0.7 km, and we report no detections in our dataset. We set a new 95% 
confidence level upper limit for the surface density Hpi(d) of KBOs larger than 1 km: Xjv(d > 



1 km) < 2.0 x 10 8 deg~ 2 , and for KBOs larger than 0.7 km T, N (d > 0.7 km) < 4.8 x 10 8 
Subject headings: Kuiper Belt, solar system: formation 



deg" 



1. INTRODUCTION 

The size distribution of objects in the Kuiper Belt is 
believed to be shaped by competitive processes of col- 
lisional agglomeration and disruption. The details of 
the structure of the Kuiper Belt size distribution can re- 
veal information on the internal structure of the Kuiper 
Belt Objects (KBOs), the history of plan et migration 
(|Kenvon fc Bromlevi[200l I Pan fc Sarill2005D. and the gas 
history in the Solar System (jKenvon fc Bromlevll2009f ). 
Large objects in the Kuiper Belt (diameter d > 30 km) 
can be observed directly in reflected sunlight. The lu- 
minosity distribution for objects larger than 100 km 
is well described by a single power-law cumulative lu- 
minosity distribution Sjv(< R) = 10 a ( R ~ Ro \ where 
£jv(< R) is the number of KBOs brighter than mag- 
nitude R per degree in the sky on the ecliptic plane , 
with an index a ~ 0.7 a nd R ~ 23 (|Fraser et al.ll2008t 
iFuentes fc Holm an 2008). This, under the assumption 
of 4% constant albedo, translates into a power-law size 
distribution n(d) oc d~ q with power index q ~ 4.5. For 
these objects the size distribution reflects the history of 
agglomeration. 

There is strong evidence for a break in the slope of 
the distribution at fainter magnitudes (smaller KBO 
sizes). Constraints on the extrapolation of a single power 
la w to magnitude greater than R ~ 35 were placed 
by iKenvon fc W indhorst ( 2001), who invoked Olbers's 
Paradox applied to the Z odiacal Background, and by 
IStern fc McKinnonl (|2000T ). who derived a slope for the 
distribution of small KBOs of q ss 3 on the basis of the 
cratering on Triton as observed by Voyager 2. A further 
hint of a break in the size distribution of KBOs is offered 
by the better-probed size distribution of Jupiter Fam- 
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ily Comets ( JFCs, iTancredi et al.l I2006L and references 
therein). These objects are likely injected into their cur- 
rent orbits from the K uiper Belt or the scattered disk 
(jVolk fc Malhotrall2008h . and it is argued that their size 
distribution would be preserved during this process. The 
size distribution of JFCs is well represented by a shal- 
low slope: q — 2.7 in the diameter range 1—10 km. 
Future Microwave Background surveys may also allow 
the setting of constraints on the mass, distance, and 
size distribution of Outer Solar System (OSS) objects 
dBabich et al.ll2007h. 

iBernstein et al. (2004) conducted a deep Hubble Space 
Telescope survey with the Advanced Camera for Sur- 
veys which led to the detection of three KBOs of mag- 
nitude R > 26.5; an extrapolation of the bright end 
power-law would have predicted a factor of about 25 
more detections for this survey. This reveals that a 
break in the power law distribution must occur at mag- 
nitude brighter than R = 28.5 (d ~ 20 km). This 
work remains the state of the art in deep direct sur- 
veys of the OSS, with a completeness of 50% at mag - 
nitu de R = 28.5. Mor e recen tly IFuentes et all (|2008h 
and iFraser fc Kavelaari (|2009l ) detected more KBOs in 
the R < 27 region of the size spectrum and better con- 
strained the slope of the bright end of the power law and 
the location of the break, R ~ 25 (d ~ 50 — 100 km, 
IFuentes et aT]|2008l . IFraser fc Kavelaarsll2009t) . 

KBOs smaller than about 30 km in diameter still elude 
direct observations. Occultation surveys are the only 
observational method presently expected to be able to 
detect such small objects in the Kuiper Belt. These sur- 
veys monitor background stars awaiting the serendipi- 
tous alignment of KBOs with the stars. The transit of 
a KBO along the line of sight briefly modifies the ob- 
served flux of the target star. At the distance of the 
Kuiper Belt (~ 40 AU) the size of the objects of interest 
is close to the Fresnel scale for visible light: this causes 
such occultation events to be diffraction dominated phc- 
nomena. The Fresnel scale is defined as F — y/ XD/2 
(|Born fc WolJll980l : lRoques et al.lll987t ). where D is the 
distance to the occulter and A the wavelength at which 
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the occultation is observed. In our survey the band- 
pass of the observation is centered near A = 500 nm 
and, at distance D ss 40 AU, F 1.2 km. Any oc- 
cultation caused by objects in the Kuiper Belt of a few 
kilometers in diameter or smaller will exhibit prominent 
diffraction effects. A diffraction pattern, characterized 
by an alternation of bright and dark fringes centered on 
the KBO, translates into a modulated lightcurve during 
the transit of the KBO along the line of sight. A unique 
feature, showing a series of wiggles, and generally a re- 
ducti on in flux, is imprinted in t he time series of the 
star ([Rogues fc Moncuguetj 120001 : iNihei et al.1 12007L see 
Figure [J) . 
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Fig. 1. — A simulated diffraction pattern (left panel) generated 
by a spherical d = 1 km KBO occulting a magnitude 12 F0V star. 
The MMT/Megacam system bandpass (Sloan r' filter and camera 
quantum efficiency) is assumed. The size of the KBO and the size 
of the Airy ring - a measure of the cross section of the event - 
are shown for comparison. The right panel shows the diffraction 
signature of the event (assuming central crossing: impact param- 
eter b = 0) as a function of the distance to the point of closest 
approach (bottom scale). The top scale shows the time- line of the 
event assuming an observation conducted at opposition (relative 
velocity v Te \ ~ 25 km s _1 ). The occultation is sampled at 200 Hz 
(dashed line), and at 30 Hz, the effective sampling rate after taking 
PSF effects into account (solid line, see Section [3). 



The overall flux reduction is dominated by the size of 
the KBO, while the duration of the event depends upon 
the relative velocity u re i and the size of the diffraction 
pattern H. We define H as the diameter of the first Airy 
ring, which it is limited by the Fresnel scale for sub-km 
KBOs and by the size of the object for large KBOs as 
follows ([Nihei et alJl200l : 



H 



(2 V3F 



D8+, 



(1) 



where the additional d6+ term accounts for the finite an- 
gular size of the star. When observing at opposition 
the relative velocity v rc \ of an object orbiting the sun at 
40 AU is about 25 km s _1 and the typical duration of an 

occultation by sub-km KBOs is about 0.2 seconds. 

Occ ultation surveys were first proposed by iBailevI 
(119761). but o nly r ecently have results been reported. 
iChang et al.1 (|2007f l conducted a search for KBO oc- 
cultations in the archival RXTE X-ray observations of 
Scorpius-Xl. They reported a surprisingly high rate 
of occultation-like phenomena: dips in the lightcurves 
compatible with oc cultations by objec ts between 10 and 
200 m in diameter. iJones et al.l (2008) showed that most 
of the dips in the Sco-Xl lightcurves may be attributed 



to artificial effects of the response of the RXTE photo- 
multiplier after high energy events, such as strong cos- 
mic ray showers. In the 90 minutes of RXTE data an- 
alyzed only 12 of the original 58 candidates cannot be 
ruled out as artifa c ts, but are hard to confirm as event s 
pones et al.1 120081 : IChang et al.1 120071 : iLiu et al.l[200l ). 
New RXTE/PCA data of Sco X-l provided a less con- 
straining upper limit to the size distribution of KBOs 
([Liu et al.M2008h . 




Fig. 2. — Megacam focal plane (McLeod et al. 2006|). A thick 
rectangle outlines a single CCD in the 9x4 CCDs mosaic. Two 
halves of each CCD (thin rectangles) are read into two separate 
amplifiers; each amplifier generates a separate output image in our 
observational mode. The x and y axis, as they would appear in a 
resulting image, are also shown. 



Several groups have con ducted occulta t ion su rveys 
in the optical regim e. iRoques et al.l (|2006|) and 
iBickerton et al.l (|2008h independently observed narrow 
fields at 45 Hz and 40 Hz, respectively, with frame trans- 
fer cameras. Such cameras allowed them to obtain high 
signal-to-noise ratio (SNR) fast photometry on two stars 
simultaneously. Both surveys expect a very low event 
rate due to the limited number of stars and the lim- 
ited exposure, and neither survey has claimed any de- 
tection of objects in the Kuiper Belt at this time 1 . An 
upper limit for KBOs with d > 1 km was derived by 
IBickerton et al.l (|2008l) by combining the non-detection 



result of the surveys oflChang et al.l (|2007h , IRoques et al.1 
(2006), and IBickerton et all (|20oj . TAOS (Taiwanese 



American Occultation Survey) is a dedicat ed automated 
multi-telescope survey ([Lehner et al.l [20091 ) that has ob- 
served a set of fields comprising ~ 500 target stars for 
over 3 years, collecting over 150,000 star hours. TAOS 
reported no detections but placed the strongest upper 
limit to date to the surface density of small KBOs 
([Zhang et al.ll2008f l. TAOS observes at slower cadence (4 
or 5 Hz) and has a relatively low sensitivity (SNR w 40 
at magnitude R — 12). For these reasons TAOS is only 
marginally sensitive to sub-km objects (with recovery ef- 



1 IRoques et alJ^2006T report 3 possible occultations from objects 
outside of the Kuiper Belt. 
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ficiency etaos ~ 0.3% at 700 m and etaos ~ 0.03% at 
500 m). 

Th e survey we report h ere was conducted using Mega- 
cam (jMcLeod et al.ll2006l Figure [2]) at the 6.5 m MMT 
Observatory at Mount Hopkins, Arizona. The use of 
Megacam in continuous-readout mode (see Section [2J 
on a field of view of 24' x 24' allowed us to monitor over 
100 stars at 200 Hz over the course of two observational 
campaigns conducted in June 2007 and June- July 2008. 
Our survey is sensitive to occultations of OSS objects 
d ~ 700 m or larger and we report no detections in 220 
star hours. Our MMT survey is designed to be com- 
plementary to TAOS and to reach smaller size limits, 
and unlike TAOS it would allow us to estimate the size 
of a detected occulting KBO. We expect further work on 
adaptive photometry and de-trending to significantly im- 
prove our sensitivity, perhaps allowing us to detect KBOs 
as small as d > 300 m. We discuss the improvements we 
are developing on this analysis in Section [7] The prelim- 
inary analysis we present here allows us to derive upper 
limits for objects d > 700 m. 

In the next section we describe the novel observational 
mode adopted for this survey. In Section [3] we describe 
the data acquired and analyzed for this paper. Details of 
the data extraction and reduction, which required cus- 
tom packages, are addressed in the same section. Sec- 
tion |4] describes the characteristics of the noise of our 
current datasets, and our noise mitigation approach. Sec- 
tion [5] describes the detection algorithm. In Section[6]we 
derive our upper limit to the density of KBOs. We also 
compare in detail the achievements of our survey to those 
of previous surveys. We draw our conclusions and outline 
future work in Section [71 

2. FAST PHOTOMETRY WITH A LARGE TELESCOPE: 
THE CONTINUOUS-READOUT MODE 

Achieving sub-second photometric sampling is a chal- 
lenge in optical astronomy. CCD Cameras can per- 
form fast photometric observations by reading out small 
sub-images, limiting the observations to yery sm all por- 
tions of the sky (e.g., iMarsh fc Dhillonl 120061). This 
is the approach adop ted by iRoques et al.l |200j|), and 
iBickerton et al.l (|2008D . who observed two stars at one 
time. Due to the rarity of occultation events, however, 
one would want to maximize the number of targets and 
the total exposure to increase the number of detections. 
TAOS achieves sub-second photometric observation on 
up to 500 targets wi th the zipper mode readout technique 
(|Lehner et alj l2009t ) . but they sample at < 5 Hz rate. 
Our continuous-readout technique allows us to observe 
the entire field of view of the camera at 200 Hz. 

Megacam, the MMT optical imager, is a mosaic cam- 
era comprising thirty-six CCDs - each with an array of 
2048 x 4608 pixels - with a 24' x 24' field of view (Fig- 
ure [2|). The standard readout speed of each CCD is 
0.005 sec/row with 2x2 binning. For this survey, we 
operated the camera in shutterless continuous-readout 
mode; that is, we kept the shutter open while scrolling 
and reading the charges at the standard readout speed, 
tracking the sky at the sidereal rate. Each star is repre- 
sented in each row that is read out of the camera, and 
the flux from a star in a row represents a photometric 
measurement of that star sampled at 200 Hz. Stacking 
each read row into a single image each star time-series 






Fig. 3. — Conventional stare mode image (one half of a CCD) of 
one of our fields (bottom panel). A series of rows from continuous- 
readout mode (center panel) from the same CCD and field, where 
the rows are stacked together in a single image. The flux profile 
of the central row of this segment of continuous-readout data is 
plotted in the top panel. 

forms a streak along the readout axis (y— axis). A small 
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TABLE 1 
Parameters of the observed fields 
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Fig. 4. — Top: power spectrum of one of our lightcurves before 
and after de-trending the lightcurve to remove noise (see Section 
13.1.21 . Bottom: power spectrum of the occultation time-series for 
a 1 km KBO at 40 AU occulting a FOV V = 12 star. 

portion of our data is shown in Figure [3] 

In this observational mode the flux from the sky back- 
ground is added continuously as the charge is transferred 
from one end of the CCD to the other, so the sky is ex- 
posed 2304 x 0.005 = 11.52 sec for every 0.005 sec 
integration on each star image (where 2304 is the effec- 
tive number of rows in each 2x2 binned CCD). In this 
mode the photon limited SNR is typically ~ 180 for an 
r' magnitude 10 star. 

When observing multiple targets simultaneously one 
can notice that the star lightcurves are affected by com- 
mon fluc tuations, or tren ds, due for example to weather 
patterns (|Kim et al.ll2008l and references therein) . In our 
observational mode, however, additional flux variations 
are caused by wind-induced resonant oscillations of the 
telescope. While the image motion along the x axis of 
the focal plane (transverse to the readout direction) can 
be resolved (see Section lB.l.ip . the image motion parallel 
to the direction of the CCD readout induces an effective 
variation in the exposure time of a star for a given row. 
These fluctuations are common to all stars in the field 
(with possible position dependencies) and therefore, in 
principle, they are completely removable. We discuss the 
de-trending of our data in Section 13.1.21 Other sources 
of noise that affect continuous-readout mode data are 
discussed in Section [4j The typical duration of a set of 
contiguous data was 10-15 minutes (after which the data 
load on the buffer would become prohibitive). For each 
amplifier, a single FITS 2 file is created wherein all of the 
rows read out during a data run are stored as a single im- 
age. For a typical run each FITS output image contains 
100K to 130K rows, corresponding to about 150-200 Mb 
of data. 

3. DATA 

We selected observing fields within 2.8° of the eclip- 
tic plane, where the concentration of KBOs is high- 
est (|Brownl 1200 1| ) . In order to maximize the number 
of targets we selected our fields at the intersection of 
the ecliptic and galactic planes (RA ~ 19 h 00 m 00 s , 
Dec ~ — 21°00'00"). We conducted our observations 
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TABLE 2 
Data set parameters. 





Start Date 


2007 June 6 


End Date 


2007 June 10 


Exposure at SN> 25 


100.61 star-hours 


Number of lightcurves with SN> 25 


990 


Number of Photometric Measurements 


7.2 x 10 7 


Start Date 


2008 June 27 


End Date 


2008 July 1 


Exposure at SNR> 25 


118.93 star-hours 


Number of lightcurves with SNR> 25 


527 


Number of Photometric Measurements 


8.5 x 10 7 
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in June- July, when our fields were near opposition (elon- 
gation angle e = 180°) and the relative velocity of the 
KBOs is highest (IRoaues et al.| [l987t iNihei et all 120071 : 
IBickerton et all I2009D . thus maximizing the event rate 
per target star. Pointing information for our fields is 
summarized in Table [T] The RA and Dec of each ob- 
served field are listed together with the ecliptic latitude 
(A) and a maximum range of elongation angles at which 
the filed might have been observed. 

We also observed control fields. These were chosen 
on the galactic plane at a high ecliptic latitude; we ex- 
pect a negligible rate of occultations by KBOs in these 
fields. These data allow us to assess our false positive 
rate. Since we report no detections the analysis of these 
fields is not discussed further in this paper. All of our ob- 
serva tions were conducted in Sloan r' filter (jSmith et al.1 
2002). A set of about 7 hours on target fields was col- 
lected in 5 half nights in June 2007 and a similar number 
of hours was collected on control fields. A set of about 7 
hours on target fields and about 6 hours on control fields 
was collected in 7 half nights in June- July 2008. Out of 
the 2007 dataset 100.61 star hours at SNR > 25 are con- 
sidered in this paper. From the 2008 dataset we use here 
118.93 star hours. Information on our dataset is sum- 
marized in Table [2j The minimum signal-to-noise ratio 
of 25, is chose n arbitrarily: 25 is t he m inimum SNR of 
the su rveys of IRoaues et alj (|2006l ) and IBickerton et ail 
(|2008ft . 3 A SNR 25 limits our sensitivity to fluctuations 
greater than 4%. An occultation of a magnitude 12 FOV 
star by a KBO of d = 400 m diameter would produce 
a 4% effect. Our efficiency tests, however, revealed our 

3 Not e however that this S NR level is obtained here fo r 200 Hz, 
whereas Roqucs et al. (2006) and Bickcrton ct al. (2008) observed 
at 45 Hz and 40 Hz. 
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sensitivity rapidly drops below 10% for objects smaller 
than d = 700 m, due to residual non-Gaussianity in our 
time-series photometric data. We discuss this in Sec- 
tion m 

3.1. Data extraction and reduction 
3.1.1. Extraction 

Custom algorithms have been developed for the data 
extraction and reduction. For each field a preliminary 
stare mode (conventional) image is collected before each 
series of high-speed runs. At the beginning of our anal- 
ysis the stare mode ima ge is analyzed using SExtractor 
(jBertin fc Arnoutsl 1 1996( 1 to generate a catalog of bright 
sources. This catalog is used to identify the initial po- 
sition and brightness of each star in the focal pane. In 
order to analyze the continuous readout data, we first 
determine the sky background for each CCD and each 
row. To do so we calculate the mean of the flux counts 
in each row after removing the measurements that are 
three ct's or more from the mean (3<7-clipping) iteratively 
until the mean converges. This removes most of the pix- 
els in the row containing flux from resolved stars. Next, 
a subset of stars that are bright and isolated is selected 
from the stare-mode catalog and used to determine the 
^-displacement of the focal plane. The focal plane is split 
into two halves, 9x2 chips each, that are analyzed sep- 
arately. We select eight stars, two near each of the four 
corners of each half-focal plane. This allows us to char- 
acterize the global motion of the targets even in the pres- 
ence of small rotational modes or spatial dependency (see 
Section |4|) . For each star (*) , and at each time-stamp 
(t), we calculate /i*(i) and <r*(t), respectively the cen- 
troid offset from the original position and the standard 
deviation of the star image, assuming a Gaussian profile. 
Note that, for a given time-stamp, flux from different 
stars will appear on different rows due to the y-positions 
of the stars on the focal plane. A 1-D Gaussian 



F± = h exp 



(a; - n*(t)f 



h 



hi 



(2) 



(where F± is the total star flux, the flux at the peak 
and Tbg the sky) is fit for each of the eight stars to each 
row of the star-streak. Thus the x-displacement fi(t) for 
all the stars in the field at time-stamp t is estimated to 
be the weighted average of the star displacements: 



J2 u*(v*{t) - n+(t )) 
*=i 

8 ' 



(3) 



where /i*(io) is the star initial ^-position and lo+ is the 
weight used for that star. 

In order to weight our average we use the correlation 
of the entire x-displacement time-series [i* with respect 
to the rest of the star set: 



UJ 



(i j) = l ^ z (f-t))(tJ-j(t) - w) ( 4 . 



w* = ^X)w(*,j); 



»?(*) «?(*) 



(5) 



where s 2 is the variance of the displacement throughout 
the duration T of the time-series. The weigh t a;* is the 
square of the Pearson's correlation coefficient ( jRicell206~il 
pag. 406), a measure of the correlation of the displace- 
ment time-series for one star with the other seven. All 
star lightcurves in the field are then extracted by aper- 
ture photometry adjusting time-stamp by time-stamp 
the center of the aperture according to the x-motion de- 
rived in this stage, and with a fixed aperture size which 
is proportional to the average FWHM in the run. 4 

3.1.2. De-trending 

The lightcurves thus extracted show evident semi- 
periodic, quasi-sinusoidal flux variations that can be as- 
sociated with oscillatory modes of the telescope in the y 
direction. In particular, a Fourier analysis generally re- 
veals two strong modes, roughly consistent among runs, 
one with period near 0.04 seconds and the other near 0.5 
seconds. Fourier spectra for one of our lightcurves, be- 
fore and after processing it, are shown in Figure |4] (top). 
Because these fluctuations affect the whole CCD plane, 
they are common to all stars and can be removed to 
achieve greater photometric precision. We now want to 
identify and remove these trends from our lightcurves, a 
process that we call de-trending. 

The gen eral algorithm we used for de-trending is de- 
scribed in Kim et al. (2008). The method takes advan- 
tage of the correlation among lightcurves to extract and 
remove common features. Since we can identify distinct 
semi-periodic modes we de-trend high and low frequen- 
cies separately (typically v > 10Hz and v < 10Hz). 

We first smooth the lightcurves, to remove all but the 
frequencies that we want to de-trend, by applying a low- 
pass or high-pass filter. We then select a subset of N T 
template lightcurves (f T ) that show the highest correla- 
tion in the lightcurve features. N T is typically about 15. 
A master trend lightcurve t is generated as the weighted 
average of the normalized template lightcurves: 

N T 

E ° 2 (frj) fr,j(t)/(Ud) 

^ = w T ^—z (6) 

3 = 1 

where the notation (/ TJ ) denotes the mean flux of f T j{t) 
over the duration T of the lightcurve, and the weight 
(T 2 (/r,j) is the variance of the lightcurve in time; r(t) 
has mean value of unity and it represents the correlated 
fluctuations in all lightcurves. 

The main trend is physically associated with an over- 
under exposure phenomenon due to global image motion 
along the y axis, which causes the effective exposure time 
to vary (see Section [2]), therefore scaling the flux. In or- 
der to remove these common trends we divide point by 
point the flux of each original lightcurve / by the trend 
master lightcurve. To improve the de-trending effective- 
ness we allow a free multiplicative factor Af (a scaling 

4 We attempted to extract the lightcurves with both fixed aper- 
ture size and variable aperture size, using the FWHM calculated 
by Gaussian fitting as a point by point estimator of the aperture 
size. The fixed aperture extraction proved to be more reliable than 
the variable aperture extraction, which induced further noise in our 
lightcurves. 
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Fig. 5. — Image motion and PSF over time: mean of the x displacements for eight bright isolated stars, at the four corners of the half-focal 
plane for two data runs (top left and right panels). PSF width from the Gaussian fit averaged over the same set of stars (bottom left and 
right). On the left we used the same run used to generate Figure [8] The arrow points to the displacement feature marked in Figure 151 
On the right the ^-displacement and the PSF width for another run, with the first 0.5 seconds shown on the left at higher time resolution. 
Note how in the second run the ^'-displacements are less prominent (note the different y scale) but the amplitude of the variability of the 
PSF is larger. 



factor) for each lightcurve as follows: 



/d,A,(t) = /(*) 



W)~ l]Af + 1 



(7) 



fd } A f is the de-trended lightcurve. 

We optimize our de-trending by selecting Af to min- 
imize the variance of the de-trended lightcurve /d with 
respect to f c = f — f s + (f s ), which is the original 
lightcurve cleaned of the frequency to be de-trended. We 
apply a high-pass (low-pass) filter to / to obtain / s if we 
want to de-trend the low (high) frequencies. Af is then 
optimized by setting: 



d 



dA 



EC/W*) - </c» 2 = o, 



(8) 



which minimizes the second moment of the de-trended 
lightcurve with respect to f c . The optimal value of Af 
can be calculated analytically. 

We set no constraints on Af, and for all of our runs 
the optimal values of Af proved to be close to 1 (which 
is what we expect in the presence of global trends) ex- 
cept for pathological cases where the flux of the star was 
buried in noise and the raw and detrended SNR were ex- 
tremely low. These lightcurves would not pass SNR cuts 
and were never considered in any of our analysis. 

Examples of the results obtained by our de-trending 
algorithm are displayed in Figures [6] and [7] In Figure [6] 
the top two panels show lightcurves for two independent 
sources in our field, and the bottom two panels show the 
same lightcurves after de-trending. Note that the top 
star is ~ 2.5 magnitudes brighter than the other and this 
is reflected in the lower SNR of the fainter source (bottom 
panel). Figure [7] shows one of our lightcurve before (top) 
and after de-trending (bottom). The raw lightcurve is 
implanted with an occultation by a d = 1 km KBO oc- 
culting a V — 9 FOV star. The diffraction feature is com- 
pletely lost in the trends and becomes evident only after 
de-trending. In the bottom panel we show the lightcurve 
detrended without allowing for the optimization factor 
Af at the top (plotted at the top at an arbitrary offset) 
and with optimization factor Af = 1.15 for the low fre- 
quencies and Af — 1.05 for the high frequencies, shown 
at the bottom. The introduction of an optimization fac- 
tor improves the SNR of the de-trended lightcurve from 



SNR = 30.0 to SNR = 30.7. For this particular 
run improvements of up to 7% in SNR were achieved by 
optimizing the detrending. 

Note that, while we used smoothed versions of our 
lightcurves to identify the trends and to optimize the 
de-trending, we do not smooth or filter our lightcurves 
to improve the SNR, thus preserving all intrinsic fea- 
tures (including potential occultations). Figure [4] shows 
the power spectrum of one of our lightcurves before and 
after de-trending it (top) . The power spectrum of an oc- 
cultation time-series generated by a 1 km KBO occulting 
a FOV star of magnitude V = 12 is shown in the bot- 
tom panel. Our de-trending greatly reduced the power 
at all frequencies: the cumulative power for this particu- 
lar lightcurve at frequencies v < 40 Hz is suppressed by 
a factor of 40. Because the oscillations are not perfectly 
correlated among our stars (see Section |4|) some residual 
power is visible. Smoothing however would would signif- 
icantly reduce the strength of the occultation features, 
that show power at all frequencies v < 20 Hz. 

4. RESIDUAL NOISE IN THE TIME-SERIES 

With a SNR > 25 we can detect fluctuations of a few- 
percent. In an 0.005 sec exposure the flux for a magni- 
tude r' = 14 star observed by Megacam is about 10 3 e~, 
which after taking into account the contribution to the 
noise of the background should lead to a Poisson limited 
SNR of about 25. While we were able to remove a large 
portion of the noise that originally affected our data, we 
typically cannot reach the Poisson-limit. We have iden- 
tified five possible sources of noise in our data: 

• Contamination by nearby sources. Overlap of 
stars along the x axis (perpendicular to the read- 
out direction) within a chip, causes reciprocal con- 
tamination in our readout mode and some of the 
stars in our fields are therefore compromised and 
excluded from our analysis. Furthermore, oscilla- 
tions of the images along the x axis causes the rel- 
ative distance between the star-streaks to change, 
which causes occasional merging. Note that while 
these oscillations are simultaneous in time domain, 
they do not occur in the same row in the recorded 
image. In each row the star images of two objects 
that are at a different y position on the CCD plane 
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Fig. 6. — Lightcurves of two independent stars in one of our 
fields. The left-hand plots show a 0.5-second chunk of the time 
series; the following 4-seconds are shown on the right at a lower 
time resolution. The top two panels show the lightcurves before 
detrending. Common modes are visible at multiple time scales. 
The bottom two panels show the lightcurves after detrending. The 
top lightcurve is the same used in Figure [4] 
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Fig. 7. — Raw lightcurve on which the occultation signature of a 
1 km KBO occulting a magnitude V = 9 F0V star has been im- 
planted (top) and the same lightcurve after detrending (bottom). 
In the bottom panel the top lightcurve is detrended without opti- 
mization (plotted at the top at an arbitrary offset) and the bottom 
lightcurve is detrended with optimization factor At = 1.15 for 
the low frequencies and 1.05 for the high frequencies. 

will not belong to the same time-stamp, therefore 
the oscillations - while simultaneous - will show a 
y offset. This is shown in Figure [51 The merging 
of streaks causes artificially high counts. Aperture 
photometry with a fixed aperture does not address 
this issue properly and fitting photometry on indi- 
vidual streaks is a computationally expensive, in- 
efficient method which is also unstable in the pres- 
ence of multiple sources close to each other. 

• Unresolved sources. Sources that are too faint 
to be visible in our 0.005 sec exposures generate a 



diffuse background. For the data in Figure [3J the 
sky level calculated as the 3a— clipped mean of the 
row counts is 140.5 ADUs. The stare mode im- 
age sky level was 48 ADUs for a 5 sec exposure, 
which would lead to a prediction of 110 ADUs for 
our 0.005 x 2304 sec effective continuous-readout 
exposure. The discrepancy is due to the presence 
of unresolved streaks associated with faint stars 
across the field. Summing all the counts in the stare 
mode image and rescaling by the exposure time of 
each row we get a number very close to the sum 
of all counts in a row of continuous-readout data. 
This contamination introduces extra Poisson noise, 
but more importantly it introduces non-Poissonian 
noise as well, since the unresolved sources are af- 
fected by the same trends the bright stars display. 
Our data shows evidence of off-phase correlation 
that might be induced by unresolved sources. 

• Positional dependency in the motion and 
trends. While we treat all of the stars in the field 
as an ensemble that moves in a solid fashion along 
the x and y axes, the image motion might also have 
a rotational component. This would lead to posi- 
tion dependencies in the motion that are not ac- 
counted for by our aperture centering algorithm. 
We have not seen evidence of dependency on the 
distance to the center of the focal plane in either 
motion or trends, but we cannot exclude that oc- 
casional rotational modes of the telescope would 
occur. Differential image motion and flux fluctua- 
tions might also be induced by atmospheric seeing. 
Both of these effects might cause the star-streak 
to move out of the photometric aperture leading 
to artificially low counts. The aperture size must 
be chosen to be such that errors due to contami- 
nation by nearby sources and errors due to streaks 
exiting the aperture are simultaneously minimized. 
Furthermore in the presence of thin clouds, vari- 
ations in the transparency might generate trends 
that would affect different sources at different times 
as the clouds move across the image. Positional 
dependencies or variations in transparency might 
contribute to the off-phase correlation of our data. 

• Scintillation. Young's scaling law (lYqunel I1967L 
iGilliland et al.lll993l iDravins et al]|1998D . 

ai = (9) 

describes the error in flux intensity / due to the 
low-frequency component of scintillation, with a = 
{AI/ 1) and where A is the telescope aperture (in 
cm), Z the angle from zenith, h is the height of the 
turbulence layer, the scaling factor ho = 8000 m, 
and AT the integration time. Competing effects 
are in place in our survey: the large aperture mit- 
igating the noise, and the low air mass contribut- 
ing to signal degradation. Note however that this 
relation holds for integration time on scales of sec- 
onds or longer. When including the effects of high - 
frequency scintillation the dependency on the aper- 
ture is expected to be steeper: 



aj oc^ 7 / 3 (secZ) : 



Cl(h)h 2 dh, 



(10) 
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where C n is the refraction coefficien t for the tur- 
bulent layer (see lDravins et al.lll997l and references 
therein) . 

Using the above equations and repres entative data 
from La Palma (|Dravins "eTaLlfl99l we estimate 
that the noise contribution from scintillation is 
(7/ < 0.01, i.e., not the dominant source of residual 
noise. As compared to the other occultation sur- 
veys the term associated to the telescope aperture 
in the SNR variance (A -7 / 6 ) is a factor 20 lower 
than the same factor for the TAOS survey, 4.5 times 
lower than the same factor for IBickertqn et ap 
(|2008l ) and 1.4 than for the IRoques et aTT T2006) 
survey. 

• Convolution of the time series with finite 
PSF. The finite size of the PSF (typically two to 
three pixels, although it occasionally was as large 
as seven) causes consecutive measurements to be 
correlated. This effect is not a source of noise per 
se, but it changes the spectral characteristics of the 
noise. The scale of this phenomenon shows up in 
an autocorrelation analysis with high power at a 
lag of about seven pixels. This is effectively a ker- 
nel convolution of our time series that smoothes 
the signal, including possible occultation signals, 
so that while we sample the images at 200 Hz we 
would expect an occultation signature to be effec- 
tively sampled at « 30 Hz (see Section [5]). Note 
that this is close to, but slightly short of, the crit- 
ical Nyquist s ampling for occultatio ns dominated 
by diffraction (jBickerton et al.ll2009t ). 

While we achieved significant noise reduction with our 
detrending, our SNR is typically a factor of two to three 
lower than the Poisson-limit. Our noise is characterized 
by high kurtosis, which is indicative of non-Gaussianity. 
Residual low frequency fluctuations (about 100 points 
period) are still noticeable in many of our time series 
(see Figure HJ . Possible improvements are discussed in 
Section El 

5. SEARCH FOR EVENTS AND EFFICIENCY 
5.1. Detection algorithm 

The signature of an occultation, sampled at any rate 
> 20 Hz, is very distinctive: it shows several fluctua- 
tions prior to the Airy ring peak, then a deep trough 
and possibly a Poisson spot feature, followed by a sec- 
ond Airy ring rise and more fluctuations (see Figure [T]) . 
The prominence of these features depends upon the mag- 
nitude and spectral type of the background star, which 
together determine the angular size, as well as the size 
and the sphericity of the oc culter, distance to the occul- 
ter, and impact parameter (|Nihei et al.ll2007l ). 

One possible approach to detecting occultations in 
our lightcurves is to take advantage of this peculiar 
shape, for example us ing correlation of templates, as in 
iBickerton et alJ (|2008l ). Given the size of our dataset, 
however, we chose to utilize a search algorithm general 
enough to capture any fluctuation of some significance, 
but which requires less computational power. We scan 
our time series for any fluctuation lasting longer than a 
duration w, and on average greater than a threshold 8 
from the local mean, which is calculated over a window 




Fig. 8. — A stare mode image (right) read out from a single am- 
plifier and a corresponding ~ 10 second chunk of high frequency 
data (left). A few bright stars and the corresponding streaks are 
indicated by letters. White arrows in the left panel point to a 
distinctive synchronous displacement feature in the data, visible 
clearly in three of the streaks, in order to focus the reader's atten- 
tion to the non-parallelism of simultaneous features in our data, 
which is due to offset in the original y position. 

W of 300 data points surrounding w. Windows w of 11, 
21, 31, 41 and 61 points were considered, in combination 
with thresholds of 0.10, 0.15, 0.20, 0.30. We define the 
local intensity I\(i) as the ratio of the flux in the local 
window w and in the surrounding window W. If the flux 
in w is suppressed by more than our threshold 9 from the 
local mean (mean over W), 

i+w/2 

E 

= < _ 

U ' i+W/2 - y ' 

E fi/ w 

j=i+W/2 

then w is considered as a candidate. This is similar to 
the Equivalent Width algorithm, which is used in spec- 
tral an alysi s, and for rare even t searches bv lRoques et al.l 
(|2006h and IWang etai] (|2009f ). Overlapping candidates 
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Fig. 9. — Efficiency plotted against the number of unidentified 
candidates (mostly false positives). Each line represents a different 
window size w, and each point represents the value of the efficiency 
at threshold 8 = 0.08, 0.10, 0.15, 0.20, 0.30, the number of false 
positives monotonically grows with decreasing 8 (larger values of 
9 on the left). All lines (all ui values) show a plateau at different 
thresholds. 

are then removed and the center of the window w that 
displayed the largest deviation is selected as a single can- 
didate event. Note that this algorithm would in most 
cases trigger two separate events for the two halves of an 
occultation on opposite sides of the Poisson spot (Fig- 
ure [T]) . These cases are later automatically recognized 
and accounted for as a single event. Different choices 
of w and 9 will produce different detection efficiency and 
false positive rate. We select an optimized subset of com- 
binations of w and 9 to be used for our event detection. 
This optimization is described in the next section. 

5.2. Efficiency 

We test the efficiency of our search by implanting sim- 
ulated occultations in our raw lightcurves. By using our 
true dataset instead of generating synthetic data we do 
not introduce any assumptions about the nature of our 
time series. We run the implanted lightcurves through 
the same pipeline as the original lightcurves: de-trending 
them and searching for significant deviations from the 
mean flux. In order to achieve better sampling of our 
efficiency the entire dataset was implanted with one oc- 
cultation per lightcurve at each KBO size we tested: 
d = 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.3, 2.0, 3.0 km, and the effi- 
ciency was assessed for each size separately. The finite 
PSF width of the star induces correlations among consec- 
utive time stamps. Given the typical PSF size in our data 
(see Figure [5]) measurements are considered independent 
if separated by more than about seven pixels. Therefore, 
to modulate the original time scries by the occultation 
signal we multiply the star flux by a synthetic occultation 
lightcurve sampled at 30 Hz. 5 

5 Since the occultation typically suppresses the flux, multiplying 
by the occultation signal reduces the noise by a factor proportional 
to the occultation flux decrease, causing us to overly suppress the 
Poisson noise by a factor of the square root of the modulation. Fur- 
thermore, sources of noise that are not proportional to the photon 
counts (such as sky background and read noise) should remain 
constant during an occultation event, but this noise is reduced by 
a factor of the flux modulation when the event is added to the 
lightcurve in this way. However, since we expect to have a very 



For the purpose of our efficiency simulations we assume 
all objects are at 40 AU, since we expect our occultations 
to be within the Kuiper Belt. There is little difference 
in the diffraction feature between 35 and 50 AU. The 
differences in spectral power between the star types do 
not impact the occultation features as observed by our 
system, so we simulate all of our occultations assuming 
an F0V type star. The angular size of the star affects the 
shape of the occultation by smoothing the diffraction fea- 
tures. It is therefore important to properly sample the 
angular size space. We find that, given the objects in our 
fields, imposing a flat prior to the magnitude distribution 
between V — 8 and V = 11 adequately samples our 
angular size range. The flat prior slightly overestimates 
the average cross section H of the events, but this effect 
is more than compensated by the loss in efficiency due 
to the fact that, for stars with larger angular sizes, the 
occultation signal is smoothed out as the diffraction pat- 
tern is averaged over th e surface of the st ar, making the 
event harder to detect (|Nihei et al.l l2007h . Overall our 
estimate of our detection rate is conservative. 

To characterize our efficiency we implant occultations 
at random impact parameters b S [0, H(d)/2]. However, 
we first want to choose the most appropriate window 
size and threshold combinations, and for that we im- 
plant occultations by d = 1 km KBOs in the reduced 
impact parameter space b £ [0, 0.3 • H(d)]. This set of 
modified lightcurves is used to optimize our parameters 
to maximize our efficiency and minimize the number of 
false positives simultaneously. Although our generic de- 
tection approach can reach high efficiency (nearly 100% 
for 1 km KBOs at zero impact parameter), it also pro- 
duces a large number of candidates, most of which are 
expected to be false positives. The combination of w and 
9 values generated efficiencies ranging between 94% (at 
w=ll and 9 = 0.1) and (at w—61 and 9 = 0.3) and the 
number of candidates ranged between and over 1000. 
Figure [9] shows the behavior of the efficiency as a func- 
tion of number of candidates. Different window sizes are 
represented by different lines and the different thresh- 
olds are marked by the points along each line. Typically, 
after a rapid increase in efficiency with the decreasing 
threshold, the efficiency stabilizes, while the number of 
candidates keeps growing: we want to choose our param- 
eters near this point, where the efficiency is highest and 
any less stringent choice would only increase the num- 
ber of our false positives. We select combinations of w 
and 9 that yield both an efficiency > 50% and a ratio 
of efficiency to candidates < 0.5. The following are the 
accepted windows-threshold combinations: (w, 9) = (21, 
0.15), (31, 0.20), and (11, 0.25). Events found in any 
run with these selection parameters were considered as 
candidates. We reached an overall efficiency of 82% at 
d = 1 km for lightcurves implanted with synthetic oc- 
cultations at varying impact parameters between and 
H/2. 

The efficiency of our search is summarized in the top 
panel of Figure [TQl as a function of KBO size. We also 
plot the corresponding effective solid angle fie (d) , defined 

high recovery efficiency for any occultation which generates effects 
> 20%, where the underestimation would become significant, we 
do not expect this effect to impact our efficiency estimation. 
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Fig. 10. — Efficiency of our survey as a function of KBO diameter 
(top). The central panel shows the effective solid angle of our 
survey and the bottom panel the effective solid angle multiplied 
d~ 2 and d~ A . 

as: 

QeC^E^^^e^U (12) 

where H(d,9*) is the cross section of the event, which 
depends on both the diameter of the KBO and the star 
angular size as indicated by (9„; w re i is the relative ve- 
locity of the KBO, which depends on the elongation an- 
gle which is close to opposition for all of our observa- 
tions; D is the distance to the occulter (assumed to be 
D = 40 AU), T* the exposure for the star target (du- 
ration of the lightcurve), and e(d, 8+) the recovery effi- 
ciency for that diameter: e(d, 9+) = 1 if the implanted 
event was recovered, otherwise. The sum is carried 
out over all of our lightcurves with SNR > 25. £le(d) 
represents the equivalent sky coverage of our survey for 
targets at diameter d, accounting for a partial efficiency. 
The center panel of Figure [TU] shows the effective solid 
angle as a function of diameter. The bottom panel shows 
the effective solid angle multiplied by bracketing slopes 
for the size distribution: <i~ 4 and d~ 2 , and it indicates 
the survey expects to see the largest number of detections 
near d = 700 m. 

5.3. Rejection of false positives 

At this stage we have more than a thousand candi- 
dates. However, most of the false positives can be re- 
moved in an automated fashion: we reject fluctuations 
that appear simultaneously in more than one lightcurve; 
those are most likely due to local weather or atmospheric 
patterns that were not corrected in the de-trending phase 
because they only affected a subset of lightcurves, and 
can be ruled out as serendipitous occultations. We also 
reject any fluctuation that does not have the right com- 
bination of depth and width. We empirically investi- 
gate the relationship between the depth and the width 
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Fig. 11. — Phase space plot showing the regions of the flux 
decrease-duration space occupied by occultations by KBOs of di- 
ameter 0.1 to 3 km, as observed through the MMT/Megacam sys- 
tem bandpass. We simulated occultations from KBOs in the size 
regime 0.1-3.0 km, and we fit the occultations with an inverted 
top-hat function with parameters A and T. The intensity of the 
gray scale reflects the number of simulated occultations with best 
fit value A and T: white areas are void of occultations. 

of an occultation by a KBO, as it is seen by our sys- 
tem, taking advantage of our simulations. To define the 
depth and width of the events we fit synthetic occultation 
lightcurves with inverted top-hat functions with param- 
eters A (depth) and T (width). Figure [TT] shows the best 
fit values A and the T for occultations simulated in the 
diameter range d = 0.1 km to d = 3.0 km, impact param- 
eters b = to H/2 and magnitude range 8 to 11 for F0V 
stars (the same set that we used for our implantation 
with additional occultations from objects d < 0.5 km). 
The shaded region represents the area of this phase space 
where at least one occultation was best fit by parameter 
values A and T (and the intensity of the shade reflects 
the frequency of A — T best fits) . We can automatically 
reject events outside the dashed polygon as incompati- 
ble with d < 3.0 km KBO occultations. 6 We are not 
sensitive to events shallower than a 10% flux drop. 

At this point the absolute number of candidates is 
small (25). The re maining candidates are inspected 
visually (using DS9, Dove k Mandel l2003| ). and the 
lightcurves are extracted with a different photomet- 
ric method (based on IRAF). All remaining candidates 
prove to be artifacts, mostly due to photometry. No can- 
didates are left after this elimination process. 

6. UPPER LIMIT TO THE SIZE DISTRIBUTION OF KBOS 
AND SCIENTIFIC INTERPRETATION 

We now comp are our limit to the si z e distribution of 
KBOs to that of lBickerton et alj (|2008D . iBickerton et all 
(2008) derived an upper limit to the surface density of 
KBOs of diameter d > 1 km. They considered the data 
obtained by their own survey t ogeth er with the data pub - 
lished bv iRoques et all (|2006f ) and IChang et alJ (|2007h . 
assuming 100% efficiency for each survey at 1 km, and ob- 
taining a total effective coverage fie = 5.4 x 10~ 10 deg 2 . 
The cross section H used to calculate fie is set to validate 

6 Note that the duration regime over which we recover events 
extend as far out as our largest window: W = 61 points or 300 ms. 
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Fig. 12. — Upper limits to the surface density of KBOs. Left panel. The dashed li ne is the best fit to the [Bernstein ct at (2004) 
survey, extrapolated to d = 0.01km. Three data-points reported by [Bernstein ct al. (2004) are plotted (HST, the faint est data from direc t 
observations). The straight line is th e TAOS upper limit to the slope of the small size end size distr ibution: q < 4.6 i Zhang et al.) 2008). 
The result by Bickcrton ct al. (2008) is shown as an empty circle (BKW) as well as the X-ray result (Jones ct al. 2008, RXTE). The upper 
limits set by our survey at d = 1 km and d = 0.7 km are shown as filled circles. The region relevant to our limit, enclosed in the square, 
is magnified on the right hand panel. Right panel. Upper limits to the surface density of KBOs, zooming in the 0.2—6 km region of the 
size spectrum where our sur vey can place limits. Symbols and labels are the same as for the le ft panel. The lower limits for JFC precursor 
populations are also shown (Lcvison & Duncan 1997; Morbidclli 1997; Volk & Malhotra 2008). 



the 100% efficiency assumption on a survey by survey ba- 
sis. Our survey adds 7.0 x 10~ 9 deg 2 to the collective fie, 
allowing us to derive a limit ove r an order of magnitud e 
stronger than the limit set by iBickerton et aTl (2008). 
Thus we set a comprehensive 95% confidence level upper 
limit on the surface density of d > 1 km KBOs at 40 AUs 
of £jv(d > 1 km) - 2.0 x 10 8 deg -2 . 

We can also derive a new upper limit for ob- 
jects as small as 700 m, where our efficiency is 
e ( d = 700 m) ~ 10%. We can set a 95% confidence 
upper limit ofEAr(d> 0.7 km) - 4.8 x 10 8 deg -2 . These 
limits are shown in Figure along with the TAOS 
model-dependent upper limit and the limit set by the 
RXTE X-ray survey. 

6.1. Comparison with the results from the TAOS survey 

Our survey aspires to be complementary to TAOS in 
that it potentially could detect objects as small as 300 m. 
However, at this stage of our work we are unable to push 
the detection limit below the TAOS sensitivity (500 m). 
Note that the recovery efficiency for TAOS at 700 m is 
ctaos ~ 0.3%, a factor of four lower than our efficiency. 

The TAOS upper limit to the surface density of KBOs 
is presented as a model-dependent limit, under the as- 
sumption of a straight power-low behavior for the small 
end of the Kuiper Belt size distribution; it is therefore 
not trivial to relate the two results, but it is clear that the 
number of star-hours a dedicated survey can collect com- 
pensates for the loss in efficiency at the small size end, 
and TAOS is able to produce more stringent limits than 
our own. Our survey would however capture the details 
of the diffraction feature with exquisite sampling, while 
the information contained in the same occultation, as 
observed by TAOS, would be greatly reduced due to the 
slower sampling. This would allow us to set constraints 
on the size and distance of the occulter, while the size- 
distance-impact parameter space is highly degenerate in 



the TAOS data. 

6.2. The Kuiper Belt as reservoir of Jupiter Family 
Comets 

The classical Kuiper Belt, the scattered disk ob- 
jects and the plutinos have all been considered in dy- 
namical simulations as possible reservoir of JFCs (see 
IVolk fc Malhotral 120081 and references therein). The 
inclination distribution of the JFCs strongly suggests 
a disk-like progenitor population, favoring the Kuiper 
Belt over the Oort Cloud. Giant planets generate long 
term gravitational perturbations that causes weak or- 
bital chaos, whi ch explains the injection of comets to 
the JFCs region (iHolman fc Wisdom 1993s: iDuncan et all 
119951: iLevison fc Duncan! Il997t ). The efficiency of this 
process depends on the dynamical characteristics of the 
progenitor family. 

Simulations of the injection process lead to lower lim- 
its on the number of progenitors, which we can compare 
with our upper limit to the surface density of KBOs. 
iBernstein et al.l (|2004l ) discussed constraints on the pro- 
genitors of the JFCs on the basis of their HST/ ACS sur- 
vey. This survey is however only sensitive to objects 
greater than ~ 20 km in diameter, while the precursors 
of the JFCs are likely to be in the size rang e 1 — 10 km. 
This is the typical observed size of JFCs (|Lowrv et al.l 
2008) and it is likely that its progenitor population would 
consist of objects of similar size (or slightly larger) than 
the JFCs themselves. 

In Figure [12] we show the lower limits to the KBO 
populations (classical belt and plutinos) and scattered 
disk derive d from dynamical s i mulat ions. We use the 
estimate of ILevison fc Duncanl (|1997l ) for a population 
of come tary precursors e ntirely in the classical Kuiper 
B elt, of | Morbidellil (l!997|) for plutinos progenitors, and 
of lVolk fc Malhotral (120081) for a progenitor popu lation in 
the scattered disk. As in IBernstein eTall (2004) we con- 
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vert the population estimates for the Kuiper Belt popu- 
lations into a surface density by assuming for each popu- 
lation a projected sky area of 10 4 deg 2 . IVolk fc Malhotral 
(2008) provide information on the fraction of time the ob- 
jects in their simulation spend between 30 and 50 AUs 
and within 3° of the ecliptic plane, and these fractions are 
used to calculate the minimum surface density of scat- 
tered disk objects expected in the region of sky typically 
observed by occultation surveys. 

We are not presently able to exclude any of these pop- 
ulations as progenitors of the JFCs. Future occultation 
surveys, with improved sensitivity, should provide valu- 
able information on the origin of JFCs. 

7. CONCLUSIONS AND FUTURE WORK 

We have devised a new observational method which 
allows fast photometry with large telescopes with stan- 
dard CCD cameras. We are able to achieve high photo- 
metric rates (200 Hz) on tens of targets simultaneously. 
The data reduction techniques for this kind of data are 
still under development. The amplitude of our noise is 
typically larger than the Poisson noise, and it displays 
obvious deviations from normality. However, we prove 
this method is suitable for gathering a large amount of 
precision fast photometric data in few observing hours. 
We present a result that lowers the upper limit set by 
similar sub-km target occultation surveys by more than 
one order of magnitude for KBOs d > 1 km, and we can 
push the upper limit to d > 700 m. We confirm the result 
obtained by dedicated Kuiper Belt occultation surveys. 

The high speed sampling achieved with continuous 
readout mode will enable the resolution of the diffrac- 
tion features of any candidate events, which is not possi- 
ble with the TAOS project due to the lower sampling rate 
they use. This will allow us to set tight constraints on the 
physical characteristics of an occulting system, possibly 
breaking the degeneracy between impact parameter, size 
and distance for sub-km KBOs. Furthermore, continuous 
readout mode enables the simultaneous monitoring of as 
many as 100 stars, which i s a dis tinct advantage over 
the su rveys of iRoaues et alj (|2006f ) and iBickerton et al.l 
( 2008), where only two stars can be sampled at a time. 
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This observational technique has proven useful in test- 
ing telescope performance and addressing guiding issues 
and it was used at the MMT to test the drive servos. 
Furthermore this observational method is a promising 
technique for ground-based high precision photometry of 
bright sources with large telescopes as it addresses many 
issues typically enc ountered in observing bright targets 
(|Gillon et al.l |20 0ct). Saturation is avoided without re- 
sorting to defocussing, it involves no overhead due to 
readout and with a camera like Megacam, with a large 
field of view, it allows the observation of many stars at 
a time, guaranteeing the presence of a good number of 
comparison stars that can be used to achieve high preci- 
sion relative photometry. 

Further improvements in SNR might be achieved: we 
are exploring a fitting photometry package that uses the 
Expectation-Minimization algorithm, treating each row 
as a mixture of Gaussians, to better separate the con- 
tribution from different sources. A possible way to ad- 
dress the contamination due to unresolved sources is to 
subtract the contribution from known unresolved sources 
(identified from the stare-mode image, see Section [5. 1 . 1 [) 
using the trends identified in the detrending phase (Sec- 
tion I3.1.2|) . Another possibility is to de-trend the 
lightcurves recursively, while allowing a variable phase 
offset. Finally, it shall be noticed that Megacam will 
become available for observations at the Magellan Clay 
Telescope, from where our target fields, at the intersec- 
tion of the galactic and ecliptic plane, could be observed 
at a higher elevation. This would help reduce the noise 
introduced by cross contamination and differential image 
motion, as the athmospheric effects we encounter observ- 
ing at high air masses would be reduced. 
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